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A SYSTEM TO MONITOR STELLAR 
IMAGE QUALITY 

Jack L. Bufton 
Optical Systems Branch 
Advanced Development Division 

ABSTRACT 

The following report outlines the development and implementation of a sys- 
tematic approach to measurement of stellar image quality in a turbulent atmos- 
phere. The concept is that of J. V. Ramsay and H. Koebier, C.S.I.R.O., Sydney, 
Australia. Simultaneous but independent measures are made of total intensity, 
size, and movement of the stellar image. The critical parameter of image size 
is deduced from a Fourier transformation of optical spread function in the focal 
plane of a six inch diameter telescope. This results in an optical transfer func- 
tion for the atmosphere and telescope optics. Experimental hardware has been 
developed and field tested at Goddard Optical Research Facility. Data is re- 
corded in the field on magnetic tape and later processed by digital computer. In 
this project most emphasis is placed on data reduction and interpretation with 
the aim of developing a capability to evaluate particular sites for their image 
quality and to make meaningful comparisons with current theories on optical 
propagation. 
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A SYSTEM TO MONITOR STELLAR 


IMAGE QUALITY 
SECTION I 
INTRODUCTION 

A device to monitor the quality of stellar images has been selected and de- 
veloped by the Optical Systems Branch. The device and its experimental program 
have been designed to meet the foUowing objectives: 

(1) Development of techniques to monitor atmospheric effects on vertical 
propagation of optical radiation 

(2) Recording of data on the optical quality of the atmosphere at particular 
sites, such as the Goddard Optical Research Facility and Mt. Hopkins, 
Arizona 

(3) Relation of data to theoretical analysis of atmospheric effects at optical 
frequencies. 

The OTTOS Study Report of June, 1967, emphasized NASA's long standing 
need for a field facility to house and operate the Agency's growing capability in 
optical tracking and communication. As a necessary beginning, the report rec- 
ommended an active program of atmospheric research and testing for site selec- 
tion and evaluation. More recently, attention has been focussed on NASA use 
of the SAO Southwest Observatory at Mt. Hopkins, Arizona, on a cooperative 
basis. Thus, plans for the Goddard Stellar Image Monitor (SIM) and other site 
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testing instruments were directed toward the specific goal of evaluating the 
Mt. Hopkins site. Primary emphasis was placed on quantitative data and mean- 
ingful comparisons with theory. 

The process of site evaluation for optical communication purposes closely 
parallels the site evaluation task in astronomy. Both cases require the input of 
information from many areas , the specification of atmospheric propagation 
characteristics being only one of the most important. Optical characteristics of 
a site include many parameters, such as cloud cover, aerosol content, infrared 
transparency, resolution of stars, and sky background. A lack of perfection in 
these measurable parameters can be attributed to three mechanisms at work in 
the atmosphere: absorbtion, molecular and aerosol scattering, and turbulence- 
induced scattering. The first two mechanisms result in effects that can be 
minimized or eliminated by proper choice of site. In addition, they are relatively 
steady-state effects, not varying significantly within a particular locality or 
length of time. The third mechanism, turbulence, results in the most serious 
optical problems. The unstable, and as yet quite unpredictable turbulent processes, 
produce random fluctuations in refractive index along the optical path. This in 
turn distorts the phase or shape of the optical wavefront and causes variations 
in irradiance across the wavefront. These distortions vary strongly over short 
periods of time and among different locations at a particular site. 

For vertical atmospheric paths, measurements of starlight yield Information 
on optical wavefront distortions. Starlight is collected by telescope optical sys- 
tems, and the quality of the image examined. This closely parallels the detection 
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of light from spaceborne laser systems by observers on the earth, an essential 
ingredient of optical communication. Astronomers have long been aware of 
degradation of stellar image quality and classed the effects as scintillation, 
dancing, pulsation, shimmer, boil and focus shift. Together these effects con- 
stitute "seeing" through the atmosphere. More compactly, they may be grouped 
into image intensity variations and fluctuations in the size, position, and sym- 
metry of the image. Until recently, measures of "seeing” have been semi- 
quantitative at best. Often, a scale of one to five has been used to distinguish 
periods of good and bad "seeing". The evaluation of a site for optical communi- 
cation requires a more exact measure of image quality. 

In the Stellar Image Monitor Experiment, image quality is determined from 
data on image intensity, size, and position. Data on image intensity fluctuations 
or scintillation is taken simply by photoelectric detection of total light in the 
image. However, size and position fluctuations of stellar images, which bear 
most directly on their quality, are not so easily measured. Fortunately, the 
data taking and analysis procedures are simplified by use of optical transfer 
functions. This is the method employed by the SIM Experiment. 

Use of transfer functions is a powerful technique especiaUy suitable for 
optical communication. A transfer function in communication theory reveals a 
system's ability to transmit information. That is, it predicts the response for a 
given input. In optical communications, the information will be a form of en- 
coding on the optical beam. The atmospheric medium will distort the encoding 
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as will the transmission and detection process. If all processes are described 
by transfer functions, they can be simply combined to describe total system per- 
formance. The optical transfer function has the additional advantage of pre- 
dicting specific parameters such as photograpliie resolution and optical hetero- 
dyne signal-to-noise ratio. 

The transfer function of an optical system is defined as the two-dimensional 
spatial Fourier transform of the image spread function. The magnitude of the 
transfer function is a measure of the reduction in contrast suffered by each 
Fourier component of the object after transmission through the entire imaging 
system. It is a function of the transform variable, spatial frequency, which is 
expressed in cycles per unit length or cycles per radian field of view in the focal 
plane. 


T(f-) 



s(x - x' ; A ) exp(- 2n i f - x) dx 


( 1 ) 


T( f ) = optical transfer function 
S(x - x' ; \) = image spread function 

For the particular case represented by starlight (point source and incoherent 
light), the optical transfer function (OTF) is just the Fourier transform of the 


image intensity distribution. ^ 


T(f) 


£ £ 


I(x) exp (- 2n i f ■ x) dx 


(2) 


x = two dimensional vector in focal plane 
I(x) = image profile 

f = two dimensional spatial frequency vector 
The magnitude or modulus of T( f ) is the quantity measured in the SIM. In 
electronic terms, useful for handling data, modulus is equal to modulation index, 
a number between zero and one. The Fourier transform inversion formula ex- 
presses image intensity profile in terms of T( f ). 


I(x) 


□ 



T( f ) exp(27r i f * x) df 


(3) 


T( f ) = optical transfer function 

Both phase and amplitude effects in the atmospherically distorted wavefront 
contribute to T(? ) . Phase changes can be thought of as tilts or bends in the 
wavefront and lead to a blurred image. Additional diffraction and consequent 
image blurring is due to wavefront amplitude effects. These effects can be con- 
sidered a random apodization of the telescope lens. The blurring transfer func- 
tion is measured in the Stellar Image Monitor (SIM) for the optical Bystem 

Appendix A contains the assumptions and mathematical analysis feinting OTF and image intensity 
profile. 
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composed of the atmosphere and telescope optics. Laboratory test data enables 
the telescope optics transfer function to be separately determined. This func- 
tion results from diffraction at the lens aperture and from the effect of any 
aberrations present. 

Besides blurring, images are also observed to fluctuate in position or dance. 
The SIM measures motion of image center of gravity as the third ingredient of 
image quality. Movements are caused by wavefront phase changes , primarily 
tilts. A transfer function for image motion may be constructed from the data 
and combined with the blurring function to yield a combined transfer function for 
the atmospheric part of the optical system. Section IV of this document will be 
concerned with relation of the transfer functions and the three data outputs of 
the SIM to atmospheric theory and engineering parameters. 

SECTION n 

STELLAR IMAGE MONITOK DESIGN AND OPERATION 

CONCEPT 

The Goddard SIM is based on a design published by Ramsay (Reference 18) 
and a technique described by Lindberg (Reference 15). Ramsay performed some 
experimental work with starlight and later Coulman (References 2 and 3) used 
the device extensively for horizontal optical propagation studies. The SIM gen- 
erates an OTF by performing a Fourier transform on the image intensity distri- 
bution in the focal plane of a telescope. The transform is generated by Bquare 


wave chopping the image at a selected but variable spatial frequency. Operation 
at different spatial frequencies allows the entire OTF of the atmosphere plus 
telescope optics to be measured. The Goddard device differs in the way spatial 
frequency is varied, in the pre-processing electronics, and in data collection and 
reduction. At Goddard, most emphasis is placed on data reduction and interpre- 
tation. Data recorded in the field on magnetic tape is later analyzed by digital 
computer. The computer is programmed to output a virtually written statistical 
record of the recorded data. A detailed account of this important process is 
given in Section in. Figure 1 illustrates the experiment concept and Figure 2 is 
a photograph of the equipment used. 

OPTICAL SYSTEM AND DETECTOR PACKAGE 

In the SIM , starlight is collected by a six inch diameter doublet lens system 
of 48 inches focal length. Figure 3 indicates the optical system components. As 
converging starlight enters the experiment package, it encounters a partially 
silvered mirror. This surface, acting as a beamsplitter, allows 30% of the light 
to pass through and reflects 60%. Transmitted light enters a penta-prism and 
is directed through a corrector plate to focus on a ring and cross-line reticle. 
The field of view corresponding to the ring's diameter is 14 minutes of arc. 
Observation of the stellar image on the reticle allows pointing of the telescope 
to insure on axis operation. Light reflected by the beamsplitter is reflected 
again by a front surface mirror and allowed to come to focus. A microscope 
objective servos as relay optics to focus an enlarged imago of star on a rotating 
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glass disk. An eyepiece and mirror assembly mounted in a sliding tube are used 
to visually check image quality and insure proper alignment of the image on the 
glass disk. The disk is composed of alternate clear and opaque pie-shaped 
sectors as detailed in Figure 3. Each sector subtends one-half degree. The 
pattern was produced at Goddard by photographic-reduction and etching from a 
machine drawn original. The disk is belt driven from a D.C. motor-servo system 
at a constant rate of 517 rpm. Speed regulation is better than . 5 %. Light pass- 
ing through the disk is thus square-wave modulated with a time frequency of 
3.1 KHz, determined by motor speed and number of sectors on the disk. The 
spatial frequency of chopping is inversely proportional to the distance from the 
center of the disk and is determined by sector width at a given distance. Since 
a clear and opaque sector determine one cycle of chopping, the spatial frequency 
in cycles per millimeter is the reciprocal of the width of two sectors. Spatial 
frequency can be made independent of lens parameters by multiplication by sys- 
tem focal length. Cycles per mm become cycles per radian field of view or 
cycles per arc second. Starlight is brought to a focus one inch from the center 
of the disk. For unity power relay optics the spatial frequency of chopping is 
2.26 cycles per mm or .0134 cycles per arc second. Higher powers effectively 
increase the spatial frequency by enlarging the image on the disk. Below are 
listed the linear increases in spatial frequency achieved with various powers of 
microscope objective. 


Table 1 


Microscope Objective Power and Spatial Frequency 


1 

! 

; Microscope Objective Power 

Spatial Frequency 
(Cycles/Arc Second) 

1 

.0134 

2.6 

.0348 

3.5 

.0469 

4.0 

.0536 

6.0 

.0803 

10 

.134 

21 

.281 

45 

.602 

60 

.803 

98 

1.31 


Light transmitted by the disk passes through an optical highpass filter and is 
received by a photomultiplier tube. The photomultiplier tube is an RCA Model 
6199 with S-ll response. The filter is a Corning glass filter number CS-72. 

The combined spectral response is approximately 1200A centered at 5000A 
(Figure 4). 

Light from a small D.C. powered lightbulb is also focussed on the disk to 
act as a phase reference source. A thin glass plate In the converging light from 
the bulb can be tilted to adjust the position of this image on the disk. Chopped 
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light is received by a PIN 10 photodiode. A field effect transistor preamplifier 
and biasing network prepares the photodiode signal for entry into the electronic 
preprocessing system. The photomultiplier signal is developed across a one 
megohm load resistance and unity gain preamplifier stage, 

PROCESSING ELECTRONICS 

Electrical signals from the two detectors are fed into an electronic pre- 
processing system 1 (Figure 5). The main function of this device is to produce 
three voltages containing information on image intensity, modulation index, and 
motion. Burr Brown integrated circuit operational amplifiers. Model 3052, are 
the primary circuit constituent. The entire instrument is D.C. powered by either 
a regulated power supply or battery providing 24 to 28 volts at approximately 
.2 amps. Solid state voltage regulators provide the plus and minus 12 volts re- 
quired for operational amplifier biasing. In addition, a programmable D.C. to 
D.C. converter provides up to 2000 volts for the photomultiplier tube. The con- 
verter is a Venus Scientific K150B rated at 1 milliamp D.C. 

The positive going input from the photomultiplier preamplifier stage is in 
the form of an amplitude and frequency modulated square wave at the chopping 
frequency, 3.1 KHz. This signal and the photodiode signal are not perfect square 
waves since image size is not negligible compared to sector width. Only the 
fundamental sine wave component of the square wave is of interest, however. 

As shown in Figure 5, the photomultiplier signal proceeds through two inverting 

^Developed under contract to NASA— GSFC by Fairchild Hiller Corporation, Bladcnsburg, 

Maryland. 

10 


1 


D.C. amplifier stages. The first has a gain of ten and the second incorporates a 
variable gain from one of twenty. At this point, the amplified signal proceeds 
two ways. The first path is through a low pass filter of the multiple feedback 
design having a corner frequency of 200 Hz. The filter's output is just the meas- 
urement of total image intensity and its scintillation. This D.C. or low frequency 
information output 1, is also used in the determination of modulation index. The 
second path is through a multiple feedback bandpass filter centered at the chop- 
ping frequency. This frequency belongs to the fundamental sinusoidal component 
in the square wave. The 3 db bandwidth of the filter is plus or minus 200 Hz. 
around the center frequency. A Transmagnetics Model 480CP5 analog divider 
receives the outputs of both filters and performs the operation A.C. divided by 
D.C. The output of the divider is the normalized chopping frequency component; 
that is, the strength of the component with respect to D.C. As the image size 
grows larger due to turbulence or as spatial frequency is increased, the A.C. 
component suffers a reduction. Size changes affect only the A.C. component 
while intensity changes or scintiUation are present in both A.C. and D.C. com- 
ponents and cancel out in the division process. The divider's output is received 
by a peak detector which produces a D.C. voltage, output 2, proportional to the 
modulation index or relative strength of A.C. to D.C. Fluctuations in modulation 
index can be foUowed from D.C. to the corner frequency of 200 Hz. If the peak 
value of A.C. is equal to the D.C. signal strength, the modulation index is 1.0. 
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Motion of image center of gravity is determined by phase comparison of 
signals at the chopping frequency, 3.1 KHz, from the two detectors. As the stellar 
image moves on the chopping disk, it produces a waveform shifted in phase with 
respect to that from the phase reference source with its stationary image. The 
preamplified reference 6ignal from the photodiode is first sent through a band- 
pass filter identical to that used in processing the phototube signal. The outputs 
of both bandpass filters are sent through limiter stages to provide good clean 
square waves. A Fairchild UA702A integrated circuit phase comparator receives 
the square waves and outputs a pulse train whose duty cycle is proportional to 
the relative phase of the signals. An RC averaging network with response from 
D.C. to the comer frequency of 200 Hz. measures the duty cycle. Relative phase 
of the two signals is just the relative position of the two images on the disk along 
a direction perpendicular to a disk radius. Thus, output 3 becomes a voltage 
proportional to one-dimensional motion of stellar image center of gravity. 

The results of electronic and optical calibration tests of the preprocessing 
system are given in Section V. These tests were necessary to insure proper 
operation of the electronics and to relate the voltage outputs to parameters 
specifying image quality. 

DATA COLLECTION 

In order to preserve the output data of the electronic preprocessing system 
in a form suitable for data analysis, recordings arc made on magnetic tape. For 
this purpose, an Ampex FR-1300 portable instrumentation tape recorder has 
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been incorporated into the system. This is a seven track machine equipped in a 
standard IRIG configuration for one-half inch wide magnetic tape. The three 
outputs of the electronic pre-processing system have bandwidths from D.C. to 
200 Hz. The D.C. requirement means a frequency modulation (FM) recording 
scheme must be used. Time signals are also recorded in order to properly 
characterize and identify the data. The NASA one per second and one per minute 
hi nary tim e codes are the ones used. These can be recorded in the direct (AM) 
mode. Recording is done at a tape speed of 7.5 inches per second. The FM and 
AM record bandwidths are from D.C. to 2.5 KHz and from 50 Hz. to 38 KHz, 
respectively. A 10 KHz sine wave is also recorded in the direct mode as a 
reference signal. Tape speed errors made in data taking can be eliminated dur- 
ing automatic data processing by using the 10 KHz reference signal. Such errors 
would result from fluctuations in power line frequency, a common occurrence in 
remote locations where power is obtained from portable generators. A complete 
set of tape recorder specifications is included in Appendix C. 

Operation of the Stellar Image Monitor Experiment involves nighttime stellar 
observations at various zenith angles, azimuths, and for various times of night. 
Two operating personnel are required: one to adjust telescope pointing and op- 
tical alignment, the other to monitor the electronics and operate the tape re- 
corder. In acquiring a star, successive use is made of a wide field viewfinder 
and narrow field, high power eyepiece. This insures operation of the telescope 
on-axis. The focussing and centering microscope objective holder is then used 
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to bring the stellar image to the correct focus and position on the chopper wheel. 
This process is monitored by the sliding eyepiece. After proper alignment is 
achieved, the eyepiece is removed and starlight enters the detector package. 

Once the photomultiplier high voltage control and D.C. amplifier gain have been 
adjusted for the proper intensity output voltage, the instrument is ready for data 
recording. 

Data taking procedures are tailored closely to the demands of data reduction. 
To this end, calibration voltages, a star identification code, and background runs 
are recorded on magnetic tape along with the data, time, and reference signals 
(see Figure 5). A voice edge track is used to note experimental parameters such 
as spatial frequency, telescope aperture, photomultiplier tube high voltage, and 
amplifier gain. Data taking is divided for convenience into four distinct cate- 
gories. A tape is the largest category. It is the contents of one reel (2,500 feet) 
of magnetic tape aal may include data from several nights , the second category. 
The third level is called a star. It includes all data taken in succession with a 
particular stellar source. Included here may be from one to ten data runs , the 
final category. 

A data run begins when the tape recorder is turned on and ends with tape 
recorder off. It is characterized by a three digit data run number and a par- 
ticular spatial frequency. Nominal operation would include six data rims per 
star-, each conducted at a different spatial frequency. Below is listed the format 
of a data run. An example is given in Figure 6. 


Table 2 

Data Run Format 


Event Number 

Event 

Length (Seconds) 


START RECORDER 


1 . 

Tape leader at ground potential 

20 

2. 

Positive six volts calibration 

10 

3. 

Ground potential 

5 

4. 

Negative six volts calibration 

5 

5. 

Star number hundreds digit 

3 

! 6. 

Star number tens digit 

3 

7. 

Star number units digit 

3 

8. 

Data 

10 to 600 

9. 

Defocus data 

5 

10. 

Background data 
STOP RECORDER 

5 


The first seven events are fed automatically to the tape recorder by a commuta- 
tor arrangement. The commutator selects the desired D.C. voltage level for the 
specified length of time. Tape leader allows time for the tape recorder trans- 
port mechanism to come up to speed and stabilize. Known calibration voltages 
are recorded to provide a voltage standard in data reduction. The three star 
numbers are the hundreds, tens, and units digits that comprise the number as- 
signed to a particular star. This number may be found in Apparent Places of 
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Fundamental Stars. Each digit is a certain percentage of plus calibration volt- 
age. For example, the star number 297 would comprise voltage levels 20%, 90%, 
and 70% respectively of plus calibration. Actual star data is taken from ten 
seconds to ten minutes. Defocus data is generated when the telescope operator 
deliberately defocuses the stellar image on the chopper wheel. This is used to 
characterize the modulation index data. Background data requires the operator 
to point the telescope so that the star is just beyond the field of view. Thus, all 
the effects of stray light, sky background, and electronic offset can be recorded. 

Data from the three outputs of the SM is recorded simultaneously in each 
data run. Time and reference signals are recorded continuously without the 
calibration voltages and star code. Below is a listing of the seven tracks of the 
magnetic tape recorder, their function, and mode of recording. 


Table 3 

Tape Recorder Format 


Track Number 

Function 

Record Mode 

1 

Spare 

FM 

2 

10 KHz reference 

AM 

3 

Image motion data 

FM 

4 

Intensity data 

FM 

5 

Modulation index data 

FM 

6 

BCD time code 36 bit 

AM 

7 

BCD time code 28 bit 

AM 


The electronic equipment described above is mounted in a standard five foot 
instrumentation rack as shown in detail in figures 7 and 8. The following items 
are included. 

1. Electronic preprocessing system 

2. Calibration voltage commutator 

3. Tektronix type 422 oscilloscope 

4. Newport Labs series 200 digital panel meter 

5. Trygon HR40-3B D.C. power supply 

6. Electro-Craft E-650 motor speed control system 

7. Ampex FR-1300 instrumentation tape recorder 

The oscilloscope and digital panel meter are used to monitor signal levels during 
data recording. Time and reference signals must be fed into the rack from an 
external time code generator. Appendix C contains selected manufactures data 
on the critical electronic components. 

SECTION in 

DATA REDUCTION 

The Goddard concept of the Stellar Image Monitor Experiment places most 
emphasis on data reduction and interpretation. What previously involved many 
tedious hours of manual data reduction from chart recordings has been perfected 
to an automated process using high speed digital computers. Data recorded on 
magnetic tape in the field is transported to the GSFC computing facility where 
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automatic routines digitize the information and read it into a computer. There 
the main program takes over, operates on the data, and outputs virtually written 
statistical records. Output is presented in tabular and graphical form for ease 
of data interpretation. 

The three outputs of SIM are randomly fluctuating voltage levels representing 
optical fluctuations in the stellar image being monitored. Similar to the atmos- 
pheric mechanisms reponsible, these optical fluctuations have a large dynamic 
range and periods as short as several milliseconds. In Figure 9 sample chart 
recordings of each of the three data tracks are presented. Clearly, one must 
turn to statistics in order to characterize the data. Consideration must be given 
to averages, variances, correlations, and power spectrums. 

COMPUTER STATISTICAL ANALYSIS 

Once data has been recorded on magnetic tape, the first step in data reduc- 
tion converts the information into a digital format. For this, the facilities of the 
Information Processing Division at GSFC are used. Analog data tapes are played 
back through an analog to digital (A/D) conversion unit and digital tapes are pro- 
duced. A/D sampling is at a rate of 4KHz. A sample is taken every quarter 
millisecond. The data tracks are sampled sequentially with the net result that 
each of the three outputs of the SIM plus a spare fourth track are sampled every 
millisecond. Fundamental sampling theorems require that the sampling rale be 
twice the highest frequency component in the data. This insures correct rep- 
resentation in digital form of the original data. Sampling 1000 times per second 


allows the data to contain frequencies as high as 500 Hz. This is more than 
enough to contain the 200 Hz bandwidths of the SIM electronics, and from pre- 
vious experience, the expected bandwidths of optical fluctuations. NASA one per 
second time code is also digitized and used primarily as a bookkeeping device. 

The digital tape contains information in binary coded decimal format and is 
ready for computer processing. 

Processing is done on an IBM system 360/91 computer through programs 
written in the Fortran IV language. Under contract to NASA-GSFC, Wolf Re- 
search and Development Corporation, Bladensburg, Maryland, is responsible 
for preparation, debugging, and operation of programs to analyze the data. The 
following description of data analysis concentrates not on specific program in- 
structions but rather on the method of statistical analysis. This method, used over 
the past two years to analyze horizontal laser scintillation and GEOS-n laser 
detection scintillation data, was first developed using a small general purpose 
computer in the Optical Systems Branch. The method is best examined by fol- 
lowing a sample data run. 

The first information input to the computer for each data track is a series 
of samples of the positive, zero, and negative calibration voltages Vj , v 2 , and 
v 3 . As mentioned previously, Vj = +6.0 volts, v 2 = 0.0 volts, and v 3 = -6.0 
volts. The computer is programmed to determine average values x t , x 2 , and 
x 3 of the samples for each calibration voltage. A linear relation is assumed to 
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exist between actual voltage v and sensed value x. 

v = A(x B) . 

Two of the resulting set of three simultaneous equations can be solved for scale 
(A) and tare (B) factors. This calibration is used later in the data run to assign 
true voltage readings to the data. It serves to eliminate inevitable gain changes 
and D.C. offsets present in the recording and reproducing process from magnetic 
tape. The computer is next presented with samples of the three digit star num- 
ber. Once again, average values are determined as well as the percentage this 
average is of positive calibration. This serves to identify the star number. The 
computer has stored a list of over 100 stars according to their star number, 
right ascension, and declination. Use of the following information allows cal- 
culation of stellar zenith angle and azimuth: 

1. Time of data run 

2. Sidereal time at Greenwich, England, for zero hours Universal time on 
the day of observation 

3. Stellar coordinates; right ascension and declination 

4. Longitude and latitude of observation site. 

Data samples then begin to enter the computer. If 20 seconds of data is to 
be analyzed, the computer will be presented with a series of 20,000 data points 
for each of tho three tracks. All statistics subsequently calculated from the data 
points will bo time averaged over the length of the data run (20 seconds). It is 
necessary to average data over a time long enough to assure good statistical 
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accuracy. Estimates based on equivalent bandwidth and resulting number of in- 
dependent samples suggest at least a ten second averaging time for SIM data. 

On the other hand, if averaging time is too long, atmospheric or other experi- 
mental conditions can change significantly, destroying the meaning of the data. 
Twenty seconds seems to be a reasonable compromise. 

Each data point is assigned a voltage level, v, according to scale and tare 
factors determined by calibration. Method of analysis then splits into two main 
directions. In the first direction, the voltage range that includes all the data is 
divided into a large number, K , of class intervals. The width, Av, of the class 
interval, i , may be chosen at will but has as its effective lower limit voltage 
resolution of A/D equipment. The following are typical parameter values: 
voltage range: , 0.0 to 10.0 volts 
K = 500 
Av = .02 volts 

The computer tabulates the number n (v . ) , of occurrences of the data record in 
each class interval and also the percentage occurrences , p (v ^ . This in effect 
establishes a probability density histogram, P (v . ) , representing in digital form 
the probability density function, P(v). 

n ( v i) 

p( v J = R - ' 

N = total number of data samples 
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1 . 


The cumulative density histogram is related to P(v ) by: 


C K) 



Av . 


This represents the total probability for voltage less than or equal to the given 
voltage, v ( . Probability density and cumulative density are related theoretically 


by: 


P(v) 


dC(v) 


C(v) 


c v 


P( v) dv . 


The computed probability density function is used to determine means, variances, 
and higher moments of the data record. 


Mean ? v 


N 





j th moment about the mean 


N 



j = 2 is the variance, c 2 . The sample mean and sample moments thus computed 
are not unique for a given sequence of N data points, since they are determined 
by choice of voltage interval endpoints and by number of class intervals. 

The parameters coefficient of kurtosis, k, and skewness, s, help to charac- 
terize the probability density function. 


m 3 

S ' (m 2 Y /2 

Coefficient of kurtosis relates to the width of the distribution and skewness re- 
lates to its symmetry. For a gaussian distribution, both k and s will be zero. 
Random voltage statistics can best be summarized in one parameter, coeffic- 
ient of variance. 


Coefficient of variance 



(4) 


This parameter has the form of a percentage modulation. It is the ratio of stand- 
ard deviation of the fluctuating signal to its average value. 

The second main direction of analysis is based on data time series; the suc- 
cession of data samples maintained in their original positions as a function of 
time. The basis tool in this analysis is the sample autocorrelation function, RC^). 


R( t ) 


1 

T-r 



v(t) v(t + r ) dt 


(5) 
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T = lag time 

T = length of data sample 
v(t) = voltage level at time t . 

The autocorrelation function for random data describes the general dependence 
of the values of the data at one time on the values at another time. An estimate 
of the autocorrelation between values of v( t ) at times t and t + r is just the 
product of the two values averaged over the length of data sample, t is the lag 
time or separation time. For computer analysis, the discrete representation is 
used. 


N-T 

R(t) = R(tAt) r N _ T v.' 

i = 0 

-t = lag number = 0, 1, 2, — m 
m = maximum lag number 
N = total number of samples 
At = interval between samples 

v ' = Vj - v = voltage level for the i ,h sample minus the average voltage 
level. 

For a good estimate of R(t), m should be chosen one tenth of N or less. 
TYPICAL PARAMETERS 
N = 20,000 
m ~ 200 

At = 1 millisecond 


Although R( t ) is a useful quantity by itself, it is more useful when applied to 
an estimate of spectral density, G( f ). 


CT 2 ( f ) 

G(f) = lim-^- 
Af-o 

Spectral density can be thought of as voltage variance per unit bandwidth. Pro- 
ceeding from the time domain, R(t) , to the frequency domain, G( f ) , requires the 
Fourier transform of R( 7 ) with respect to t. 


•CO 

G(f) - 4 R(t) cos (2v fr) dr 


(6) 


In discrete form for computer analysis, this becomes: 

V -1 v-tf 

G(f) = 4At > R(-fAt) cos “f— 

-t-o 

f = frequency (Hz) 

f c = 500 Hz = frequency cutoff of data 
The equivalent bandwidth or frequency resolution in this analysis is B s . 

= IT 

For the previously assumed parameters B c = 5 Hz. The raw estimate spectral 
density requires smoothing because the variability of the estimate does not de- 
crease with increased sample size or record length. Smoothing or weighting 
nonuniformly can be accomplished by multiplication by the H anning lag window. 
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MODULATION INDEX DATA ANALYSIS 


The modulation index track results in a voltage proportional to the modulus 
of the optical transfer function for the particular spatial frequency used in data 
taking. Voltage statistics are appropriate and require only a conversion factor 
to modulus statistics. Laboratory tests generate a plot of modulus versus modu- 
lation index voltage. Modulus is just the depth of modulation of the square wave 
input to simulate the waveform from the detection of chopped starlight. One 
hundred percent depth of modulation corresponds to a modulus of 1.0. As the depth 
of modulation is varied, the analog divider and peal: detector respond with a 
change in voltage. Modulus values obtained from this curve are corrected for 
the effect of electronic and optical background. This is done with average 
voltage and average background data obtained from the image intensity 
track. The average value of defocus voltage serves as the zero reference for 
modulation index voltage. The variance of modulus is read directly from the 
calibration curve. Data reduction is incomplete without the time series informa- 
tion provided by autocorrelation function and spectral density. This analysis can 
be performed directly on modulation index voltage. Secant of stellar zenith angle, 
azimuth and spatial frequency serve to characterize the data run. 

IMAGE MOTION DATA ANALYSIS 

The image motion track may be handled entirely in terms of voltage statistics, 
since optical phase, image motion, and output voltage of the electrical phase 
comparator enjoy a linear relationship. Image motion in the telescope focal 


plane is the result of random deviations or tilts in the angle of arrival of various 
parts or all of the incoming wavefront. These tilts are actually optical phase 
differences following a Gaussian distribution. Stellar image motion with respect 
to reference light image causes a change in the relative electrical phase 
of waveforms from the two detectors. Suitable laboratory calibration yields a 
conversion factor in volts per arc second for each spatial frequency. Voltage 
statistics are thus converted to image motion statistics. The conversion to 
average image position is meaningless in data reduction, since it is not important 
to know precisely where the image is but how it fluctuates with respect to a fixed 
point. The important parameter is image motion variance, c/. Either this 
number or rms value, ■' ,, , is plotted versus secant zenith angle of the star. 
Numbers on image motion must be multiplied by /2 to give total two-dimensional 
image motion. The measured motion is only the component perpendicular to the 
chopper disk radius. This analysis assumes isotropic image motion. Spectral 
densities of image motion can be scaled directly from voltage spectral densities. 

At the conclusion of computer analysis of stellar data for a particular night, 
a data summary sheet is prepared. The following information is tabulated for 
each data run. In addition, data for different stars is combined to yield plots of 
log e (CIV) versus log, [secant (zenith angle)] and phase variance versus secant 
(zenith angle). A modulation transfer function is plotted for each star observed 
at four or. more spatial frequencies. This places the data in a form suitable for 
interpretation. 
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Table 4 

General Parameters for Each Data Hun 


1. Data run number 

2. Star name and number 

3. Zenith angle 

4. Secant (zenith angle) 

5. Log e (secant (zenith angle)] 

6. Azimuth 

7. Time interval for data run 

8. Spatial frequency 

9. Aperture diameter 

Table 5 

Statistical Summary for Each Data Run 
INTENSITY TRACK 

1. CIV = cr y j\ - coefficient of irradiance variance 

2. of = log amplitude variance 

3. v = signal mean 

4. log 5 (CIV) 

5. v B = background mean 

MODULATION INDEX TRACK 

1. o v/ V 

2. v = signal mean 

3. v D = defocus mean 

IMAGE MOTION TRACK 

1. v = signal moan 

2. o ra 2 = motion variance 

3. = motion rms 


SECTION IV 

DATA INTERPRETATION 


BACKGROUND 

Data interpretation is extremely important for the SIM experiment. Data 
points are largely dependent on many experimental parameters such as tele- 
scope aperture size, local weather, length of data run, wavelength of light, and 
stellar zenith angle. It would be a hopeless task to control all parameters in any 
program of experiments. Thus data taken in this experiment must be examined 
in light of recorded parameters, theoretical analysis, and experimental work of 
other researchers. Such an examination must be directed toward useful informa- 
tion. Engineering numbers must be arrived at which serve to characterize the 
data and predict the performance of other optical systems in the presence of 
atmospheric turbulence. In particular numbers arc required to specify the in- 
tensity, blurring, and motion of an image formed in the optical system. 

In recent years much theoretical and experimental work has been done on 
the effect of atmospheric turbulence on optical propagation. To date there is no 
unified theory applicable to this problem , although several treatments have been 
partially successful in predicting experimental results. The usual starting point 
is an approximate solution to the scalar wave equation, although some researchers 
have appealed directly to two point correlation statistics. All approaches attempt 
to arrive at expressions for log amplitude and phase of the optical wavefront in 
a plane perpendicular to the direction of propagation through the turbulent 
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square difference of refractive index. 


atmosphere. Since amplitude fluctuations are multiplicative rather than additive, 
log amplitude is modulated in an additive fashion as is phase. As a consequence 
of the central limit theorem log amplitude and phase follow a normal or Gaussian 
distribution. Statistics of log amplitude and phase are means, variances, and 
structure functions. The latter is a particularly useful way of relating the in- 
formation. It is the mean square difference of values of the physical quantity at 
two points separated by a distance, r. Log-amplitude, i , and phase, f., structure 
functions are defined as follows. 


Df, (r) = -^W] 2 } 

(8) 

mo = ([*('.) 

(9) 


Together they form the wave structure function. 

D„(r) = D^(r) + D 0 (r) (10) 

Theory relates optical statistics (log amplitude and phase statistics) to the 
statistics of the intervening atmosphere. For this a certain knowledge of the 
physics of the turbulent process is needed. Statistics of turbulent mixing of air 
parcels are given in Reference 14 by the Kolmogorov theory. Different parcels 
of air are at different temperatures and presumably in a certain subrange of 
sizes, all energy dissipation is through turbulent mixing. Parcel sizes range 
from an inner scale, , on the order of millimeters at the earth's surface, to an 
ill-defined outer scale, tens to hundreds of meters. The important result of this 
theory is a refractive index structure function, D N (r; h), expressing the mean 
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D N (r;h) = C N 2 (h)r 2/3 (11) 

C N 2 ( h ) , the refractive index structure constant, is a function of altitude above 
the earth’s surface as well as atmospheric conditions. It is essentially a meas- 
ure of the strength of turbulence. The turbulent profile, C N 2 (h) , versus altitude, 
appears in all equations derived from the theory. Data provided by Hufnagel 
(Reference 8) gives the most accurate results to date on average turbulent pro- 
file. This data was compiled through a combination of stellar scintillation and 
meteorological measurements. As seen in Figure 10 the profile decreases quite 
rapidly from a maximum near ground level. The primary distinguishing feature 
is a sharp peak near 12 kilometers. This has been attributed to the region of 
high turbulence at the tropopause. 

Rather than proceed with a detailed treatment of solutions to the wave equa- 
tion and their relation to turbulence spectra, end results will be quoted directly. 
Only these results which apply to ground based observation of starlight will be 
presented. Expressions will be arrived at which predict the data to be obtained 
from SIM measurements of stellar irradianee fluctuation, image size, and image 
motion. 

IRRADIANCE STATISTICS 

Irradianee fluctuations have been treated by Tartarski (Reference 19). He 
has arrived at an expression for log amplitude variance, ’ i 2 , of starlight. This 
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is just four times log irradiance variance. 

C N 2 (h)h s Slh (12) 

h = altitude 
9 = zenith angle 
k a 2n/\ - wave number 

C N 2 (h) = refractive index structure constant 
This expression reveals the dependence of scintillation of wavelength, zenith 
angle, and strength of turbulence. Wavelength dependence is weak and in the 
visible spectrum amounts to only a small correction factor. In the SIM the 
optical bandpass is approximately .1 microns. For purposes of scintillation this 
can be considered monochromatic. On the other hand zenith angle dependence 
is quite strong. It expresses the relation of scintillation to the amount of air 
starlight must traverse in order to reach the telescope. Air mass is proportional 
to secant (zenith angle). Dependence on turbulence is expressed by the integral 
over turbulent profile, C N 2 (h) . 

Log amplitude variance is just the statistical parameter computed from SIM 
intensity output. It is the variance of natural log of the ratio of instantaneous 
intensity in the stellar image to mean intensity. The other useful statistical 
parameter is CIV, coefficient of irradiance variance. A straightforward 


. S6k 7/6 ( sec 9 ) 11/6 
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statistical analysis relates CIV to log amplitude variance. 


irradiance I - (i) exp(2f) 

log amplitude -C = "" 2 ~ log e 1 /( 1 ) 

irradiance variance " o * (('-«)’> 

log amplitude variance - (('-<->)■) 

c, 2 <|<l)exp(2-f)- (<I>ex P (2^))]^) 

' (I) 2 (exp(2',)))^> 


Since { is a gaussian random variable 


<exp( 2f )) •- exp l(') ( 2) 2 ((’ “ < 2 >) ’)] 

exp 2(f) + 2^-f - (•?)) 

- exp^2(f) 4 2o/ 2 J 

°i - - ( { ) 

(exp( 2f )) ' I 

u \ - (I) 2 <[exp(2t)-l] 2 ) 

= (l) 2 (exp(4()- 2 expf 2-f. ) * l) 

(exp(2f)) 1 

(exp( 4f )) = exp [ 4 (f) 4 8.. /j 

= exp j4o 7 j 
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This effect is attributed to averaging over fluctuations from different areas of 
the aperture. Both Fried (Reference 4) and Tatarski (Reference 19) predict 
aperture averaging, which not only reduces scintillation but changes its zenith 
angle dependence. In this analysis the correlation distance for irradiance 
fluctuations becomes important. 

correlation distance - H 0 sec 6 
8 = zenith angle 
\ - wavelength 

H 0 = scale height of atmosphere 
= approximately 15 km. 

If telescope diameter, D, is- greater than the telescope aperture averaging. will 
take place: while the region D/Va. see 3 < l corresponds to little or no 
averaging. The effect of exceeding the correlation distance has been measured 
experimentally by Protheroe (Reference 17). He obtained a 2.5 times reduction 
in CTV in going from a one inch diameter telescope to a six inch diameter. This 
corresponds to a reduction in a/ of 6.25. For. the six inch diameter telescope, 
he observed the following dependence on zenith angle. 

af a sec 3-4 6 
CIV “ sec 1 ' 3 8 

The calculated correlation distance was around three inches. Protheroe's actual 
CIV values must be corrected for his experimental technique in order to arrive 
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at CIV for zenith viewing. When this is done the following average values can be 
presented, which in turn predict average * values. 


Table G 

Average CIV and oj? Values 



erv 

2 

i 

Summer 

.14 

.0049 

Winter 

.16 

| .0004 


Returning to equation 12 it is noted that the integral over C N 2 (h) includes 
a factor h s 6 . Lower altitudes where C N 2 (h) is the largest do not have an im- 
portant effect on c^ 2 . It is the altitudes where h s 6 C N 2 (h) takes on its larger 
values that figure most prominently in scintillation. There are good practical 
reasons why scintillation should originate at high altitudes. 

(1) The most likely mechanism for generating temperature fluctuations 
(and thus refractive index fluctuations) is through the combined agencies 
of wind shear and non-adiabatic lapse rate present at the tropopause, 

(8 to 14 km). (Reference 13) 

(2) Consideration of turbulent element sizes and effective scattering angles 
show that altitudes of 8 to 14 km are required for the interference pat- 
terns known as scintillation to be produced on the ground. 

The frequency spectrum of irradiance is determined by the same factors of 
correlation distance, wavelength, zenith angle, and aperture size. In addition the 
velocity spectrum of upper altitude winds becomes very important. Increases in 


either wavelength, zenith angle, or aperture size are reflected in a decrease in 
magnitude of scintillation high frequency components. These effects would be 
observed in spectral densities of stellar image intensity and log amplitude. 
Tatarski (Reference 19) has developed a theoretical analysis of spectral density 
parameter dependence. 

The physical basis of spectral density aperture dependence can be under- 
stood by considering the structure of irradiance in the plane of the telescope 
aperture. The structure consists of a continuous size spectrum of bright and 
dark patches changing with respect to time and sweeping across the aperture. As 
the aperture becomes much larger than the correlation distance the small patches, 
which have dimensions on the order of correlation distance and are responsible 
for high frequency components, produce less of an effect due to aperture averag- 
ing. Larger patches and therefore lower frequency components do not suffer as 
much averaging. The rate at which the pattern is swept across the aperture is 
determined by upper altitude wind profile. Wind velocity is strongly peaked in 
the 8 to 14 km region of the tropopause, previously asserted to be the region of 
scintillation origin. This comes as no surprise since turbulence is readily as- 
sociated with high flow velocities. Typical velocities at the tropopause are from 
20 to 100 meters per second and the flow is often referred to as a jet stream. 

On the basis of the above physical considerations a proportionality to irradiance 
frequency spectrum can be asserted. 

f a v n/ ^\ H 0 sec 0 
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\ H 0 sec d = correlation distance 


v n = wind velocity normal to path of propagation 
This relation can be used to define a characteristic frequency, f ,, , equal to the 
right hand side above, f c predicts the maximum frequency of scintillation to be 
observed if v n is the maximum wind velocity, f c can be measured experimentally 
and then compared with wind speed and zenith angle data. 

OPTICAL TRANSFER FUNCTIONS 

To analyze atmospheric effects on modulation index or the modulus of the 
optical transfer functions , it is necessary to return to equation 10 for the wave 
structure function. 

D„ ( r ) = Rf. ( r ) + ( r ) 

Tatarski has related the refractive index structure function, D N ( r; h), from the 
Kolmogorov theory to the wave structure function and arrived at the following 
functional dependence. 

D„(r) = 2.91 k 2 r 5 ' 3 |c N J (z) dz 

•'path of props gat ion 

For stellar observations the path of propagation, z , should be expressed in 
terms of altitude, h, and zenith angle, 



D„(r) - 2.91k 2 r 5 3 seed C N 2 (h) dh (15) 

k = 2” - K = wave number 
r = \f 

f = spatial frequency 
r = zenith angle 
h = altitude 


Different theoretical approaches by Fried (Reference 6) and Hufnagel and 
Stanley (Reference 9) have resulted in the same expression for modulus of OTF 
in terms of the wave structure function, 



The angular brackets represent time averages short with respect to the effective 
period of atmospheric fluctuations. M o f ( f ) , the modulus measured in SIM, is 
in reality averaged over a time equal' to the reciprocal bandwidth of the modula- 
tion index electronics (sec Appendix B). The bandwidth Is 200 Hz; thus the 
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averaging time is approximately five milliseconds. The raw modulation index 
data from the SIM electronics is thus a measure of fluctuations in short term 
blur circle. Computer generated statistics are average, variance, and higher 
moments of the short term blur and do not include effects of image motion less 
than 200 Hz. Blur data in the form of a transfer function contains contributions 
of the telescope optics and the atmosphere. Laboratory test data provides a 
telescope optics transfer function. This is divided into the measured M o % ( f ) to 
arrive at the atmospheric contribution. 


j-D w (0 

1 

f* 

Ai 

,} 

l 

L 




m l (f) 


< M A< f )> sh „ r[ (17) 


M L (, f) = telescope optics transfer function 
For short term blur the atmospheric OTF is not entirely independent of lens 
parameters. Equation 17 reveals a dependence on lens diameter. However, the 
effect of lens aberrations are removed by the division of transfer functions. 

Functional dependence in Equation 16 involves wavelength, zenith angle, and 
strength of turbulence, although the specific form is different than that for ir- 
radiance statistics. Modulus dependence on each parameter can be verified ex- 
perimentally by holding constant all but the parameter under test. The strength 
of turbulence factor now is an integral over only the refractive index structure 
constant. For irradiance the integral was the 5/ 6 ths moment of C N 2 ( h ) with 
altitude. Comparison of measured results of the two integrals can result in a 
crude knowledge of the vertical profile of turbulence; i.e., the dependence of 
C N 2 (h) on h. 


The atmospheric modulus can be written, exp j- cf 5/3 [l - t (\f/D) 1/3 ]J 
where f is spatial frequency in cycles per arc second and C is a factor contain- 
ing the parameter dependence. This predicts the shape of atmospheric OTF. 
For spatial frequencies much less than the resolution cutoff, DA the shape is 
controlled by the factor, exp (- cf s/3 ) . As the resolution cutoff is approached 
the factor exp [c e (\f 5 /T>) 1/3 J becomes dominant. Using the expression 
exp (- cf s ' 3 ) the dependence of image profile on the parameters contained in C 
can be derived. The two-dimensional Fourier transform result derived in Ap- 
pendix A provides the necessary relation. 


I(s) 


r- 


M(<0 J 0 (2* so) 


= spatial frequency (cycles/arc second) 
s = arc seconds 
Now, substitute in: 


M(oj) - exp (- c<w 5/3 ) 

I(s) “ 2 I to exp (- cc,; 5/3 ) J 0 ( 2vso ) t\f, 

Jo 

^ 2 J a> exp (- coj 2 ) J 0 (2r/s^)d&. 

Jo 


exp l- 


c A 2 
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This is a gaussian shape with variance o s 2 . 


IMAGE MOTION 


c 

2tt 2 


or 


f T s 2 = 5.8 1/3 sec 9 I C N 2 (h) dh 


1 c " 2 


2.4 A -1 ' 6 ( sec 9 ) 1/2 


C N 2 (h) dh 


a , is the rms value of I(s) , i.e. the rms of image size. It is seen to be propor- 
tional to the negative sixth root of wavelength, the square root of zenith angle 
and the square root of integrated turbulence. 

Rather than image profile a specification of optical resolution, R, through 
the atmosphere can be made. 


Fluctuations in the position of image center of gravity result from atmos- 
pherically induced changes in the phase or shape of the incoming wavefront. 

Most changes (90%) can be attributed to random linear phase shifts or tilts in 
the wavefront. Tilts cause image displacement. The other principal wavefront 
degradation Is quadratic in form and leads to focusing errors; thus contributing 
to blurring of the image. Tatarski (Reference 19) has related wavefront tilts to 
phase structure function by examining a two point correlation or interferometer 
experiment (Figure 11). 

Consider two point detectors A and B separated by a distance d. Incoming 
light represented by a wavefront tilted with respect to d is brought to a focus 
after passing through the point detectors. The angle of tilt, Aa, radians, is equal 
to the angular displacement of image in the focal plane. The phase shift, A <j>, 
between A and B can be expressed in terms of tilt angle. 


R 


J. M L< f > 


(18) 


Resolution is just the area under the OTF curve. It is equivalent to the Strehl 
definition of image quality. 


A 0 = ~ - kd Aa 

The mean square difference of phase is just the definition of phase structure 
function. 

<(Ad0 3 ) = ((^b-^a) 2 ) = tVW 

- k 2 d 2 ((Aa)^ 

■■■ <( Aa)2 > = -j^r as) 
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is a starting point for the usual tests of parameter dependence. Strangely enough 
Equation 21 predicts that image motion is independent of wavelength and has 
only a weak aperture dependence D' 1/3 , for apertures greater than the inner 
scale of turbulence (a few millimeters). The table below lists aperture dependence 
for various diameters. 

Table 7 

Image Motion Aperture Dependence 


Aperture Diameter 


D -l/3 

(inches) 

(cm) 

(cm' 1 ' 3 ) 

1 

2.54 

.787 

3 

7.62 

.509 

6 

15.2 

.404 

12 

30.5 

.320 

I 60 I 

152 

.188 


Based on this analysis variance of image motion present in a six inch diameter 
telescope would be reduced only by a factor of about two for a sixty inch diameter 
telescope. Image motion may be a factor of concern in even the larger telescopes. 


Image motion dependence on turbulence is expressed by the integral over 
C N 2 (h) . The integral appears in the same form as that for the modulus of the 
optical transfer function. Since C N 2 (h) appears alone in the integral, turbulence 
strength at each altitude is weighed equally. Lower altitudes where C N 2 (h) takes 
on its largest values contribute most heavily to the integral. Due to the sensitivity 
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of image motion and size data to local lower atmospheric conditions, experi- 
mental results of various researchers are difficult to compare accurately. 
Quoted below are average nighttime results reported by Hosfeld (Reference 7) 
and Johnson et al. (Reference 12). These may be used for general comparison. 


Table 8 

Average Image Size and Motion 


i 

j Rms Motion j 

Rms Size [ 


i 

(Arc Seconds) 

(Arc Seconds) 

Hosfeld 

.62 

1.3 

1 Johnson et al. 

.34 

— 


References 9 and 10 detail the analysis of image motion in terms of a 
transfer function ( f )\ . The formula arrived at: 

'long 

= exp(-47T 2 f 2 <r/) (22) 

f = spatial frequency (arc seconds) 
ct s 2 = one-dimensional image motion variance 
is a good approximation, but not rigorously correct. If short term blur and 
image motion are assumed uncorrelated, they can be simply combined to yield 
an expression for the long term blur. 

= <M A (f)>, horl <M„(f)> long (23) 

This is the OTF to be expected when the measurement is made over a long time 
or the telescope has an extremely large aperture; hence little image motion. 


Fourier inversion of (M a ( f )) long leads to long term image profile; including 
all the effects of blur and movement. This is the quantity proportional to optical 
heterodyne S/N (Reference 5). 

SECTION V 

LABORATORY AND FIELD TESTS 

ELECTRONIC TESTS 

Prior to deployment of the Stellar Image Monitor Experiment in the field, 
a series of laboratory tests and calibration runs were performed. The first 
system component to be tested was the electronic pre-processing system. 

Figures 12 and 13 are measured frequency and phase response of the bandpass 
and lowpass filter sections. Bandwidths are 3100 Hz ± 200 Hz and 200 Hz, re- 
spectively. Frequency response was measured by using a sine wave signal 
generator in place of inputs from photomultipUer and photodiode. Phase response 
was measured on an oscilloscope. Filter input and output were applied to the 
vertical deflection and horizontal sweep respectively. The resulting Ussajous 
pattern displayed on a phase reticle gave readings of phase difference. Gain and 
linearity tests of the amplifiers were also performed. This served to calibrate 
the image intensity track. 

Modulation index caUbration required a measure of combined response of 
analog divider and peak detector. Sine waves at the chopping frequency were 
input in place of signals from the photomultiplier. Modulation index (peak A.C. 
divided by D.C.) of this input was measured at the input terminals of the analog 
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divider. Measured values of peak detector output voltage for the given modulation 
index or % modulation are shown in Figure 14. This curve applies for D.C. volt- 
age in the range -1.0 to -3.0 volts. It should be noted that linearity is not an 
important consideration in this output since the expected fluctuations will be 
small compared to their average value. Percent modulation determined from 
this curve must be multiplied by 

D.C. + BGND / D.C. GA!n \ 

D.C. \ A . C . GAIN / 

to correct for filter gains and any D.C. offset voltage. The symbol BGND, or 
background, represents this offset. D.C. + BGND information is available from 
the image intensity track. Filter gains were measured during calibration and 
their ratio is: 

D.C. GAIN 

A . C . GAIN ' 1 025 

OPTICAL TESTS 

Calibration and testing of the image motion electronics employed an artificial 
star source, Figure 15. A point source at the focus of an f/10 twenty inch diam- 
eter collimating mirror represented the small angular size and parallel light of 
a star. The point source was either white light or .6328 micron red light from 
a Spectra Physics 131 helium-neon laser. Light from an incandescent lamp 
filament focused on a variable size pinhole served as the white light point source. 
In both cases a further reduction in spot size was achieved through microscope 
objectives. Source sizes were determined as follows. 
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White Light 



Therefore each thousandth of an inch in the focal plane corresponds to one arc 
second. Source sizes were kept smaller than the theoretical angular resolution 
of six inch diameter optics. 

Collimated light produced by this test set-up was allowed to enter the six 
inch diameter optics of the Stellar Image Monitor after traversing about 40 feet 
of laboratory air space. The point sources were mounted on an optical bench 
equipped with a translation slide. Movement of the point source in increments 
of .001 inch in the focal plane could be measured there and produced the cor- 
responding number of arc seconds deviation. The D.C. voltage output of the SIM 
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phase comparator was recorded as a function of this arc second displacement. 
The output of the phase comparator was not a pure D.C. voltage but contained 
frequency components due to room vibration, air path turbulence, and chopper 
pattern centering error. The resultant plot, Figure 16, yields the desired cali- 
bration factor in volts per arc second. This factor depends linearly on spatial 
frequency. For the spatial frequency used in the calibration the rms error was 
approximately .1 arc second. It should be noted that the calibration curve is 
cyclic with a linear response in the mid range of voltages. Thus in actual data 
recording the average value of phase comparator output should be maintained at 
approximately 3.5 volts. A check on spatial frequency was made by noting the 
number of cycles of. the response curve per arc second. The results agreed 
with microscope objective power being used and theoretical spot size 
calculations. 

Using the artificial star system described above optical transfer functions 
for the SIM in white light and red light were generated (Figure 17). This was 
done by varying microscope objective power in the SIM. Measured OTF shows 
a degradation of performance in white light while monochromatic red light per- 
formance is close to the diffraction limit. 1 Chromatic aberration causing this 
effect must be present in the SIM optics since the collimating system is 
reflective. 

! Scr Appendix A foe OTF of a diffraction limited optical system with six inch diameter circular 
aperture. 


FIELD TESTS AND DATA 

All field tests to date have been conducted at the Goddard Optical Research 
Facility located approximately 4.5 miles from GSFC in an open field. The SIM 
was placed on a 12 inch high concrete slab approximately 15 feet from a 20 by 

'i 

20 foot trailer. This building supplied the necessary power and support capabil- 
ity. Time and reference signals were available from an Astrodata time code 
generator located in an adjacent building. The telescope mount was roughly 
aligned in the North-South direction. Standards alignment procedures were fol- 
lowed to set the polar axis for the latitude of 39 degrees 1 minute and to make 
fine adjustments in North-South pointing. Data runs were taken on various clear 
evenings during fall-winter 1968-1969. Of the data analyzed by computer, rep- 
resentative statistical records are presented here. j 

Figure 18 contains chart recordings of the three tracks of the SIM for data 
run 118. Experimental parameters that characterize the data run are listed be- 
side the chart recordings. Timing marks at one second intervals are located 

I 

along the top. Appropriate calibrations are indicated for the vertical scale on . 

each recording. The modulation index track' has indications of percent modula- 
tion and the image motion track is calibrated in arc seconds. The computer 

generated statistical record is presented In table 9. The three salient results i 

are: 

log amplitude variance = .018 
OTF modulus = .78 

image motion variance = .18 (Arc Seconds) 2 
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Table 9 


Figures 18 through 20 summarize data taken during the following nights. 


Statistics for Data Run 118 


GENERAL PARAMETERS: 

1 . DATA RUN 118 

2. STAR: ALDEBARAN NO. 168 

3. ZENITH ANGLE (ZA) = 29.7° 

4. SEC (ZA) = 1 .15 

5. LOG„ [SEC (ZA)] = .141 

6. AZIMUTH =112° 

7. 20 SECOND DATA RECORD: 301-08-46-36—56 

8. SPATIAL FREQUENCY = .047 CYCLES/ARC SECOND 

9. APERTURE DIAMETER = 6 INCHES 

STATISTICAL SUMMARY: 


INTENSITY TRACK 
1 . CIV = .26 

2. a? = .018 

3. v =2.84 VOLTS 

4. LOG. (CIV) = -1.34 

5 . v B = .40 VOLTS 

MODULATION INDEX TRACK 
1 . cr y /v = .05 

2. V =3.9 VOLTS 

3. v D = .47 VOLTS 

IMAGE MOTION TRACK 
1 . v = 4.2 VOLTS 

2. o,* = .18 (ARC SECONDS) 2 

3. o m = .43 ARC SECONDS 


Table 10 

Dates of Stellar Observations 

Night 


Date (Zulu) 

Local 

269 


25 SEP 68 

24 SEP 68 

292 


8 OCT 68 

7 OCT 68 

301 

, 

17 OCT 68 

16 OCT 68 

324 

| 

19 NOV 68 | 

18 NOV 68 


Included here are statistics for 58 data runs. This amounts to an analysis of 
over 110,000 data points. 

In Figure 19 irradiance statistics are presented in a plot of log e (CIV) versus 
log r [secant (zenith angle)]. All data was taken with a six inch aperture. Pro- 
theroe’s experimentally observed zenith angle dependence is plotted as a straight 
line through the data points. There is good agreement with this dependence for 
zenith angles less than 60 degrees. For larger zenith angles the effect of stellar 
azimuth becomes important and data scatter is evident. Vertical axis intercepts 
yield the values of CIV for zenith viewing; 11% and 21% respectively on the two 
nights. These values compare with Protheroe's average values of 14% and 16% 
for a six inch aperture. CIV may be converted to log amplitude variance, , 
using the approximation: 

CIV 2 * 4o^ 2 
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Spectral densities of irradiance are shown in Figures 20, 22, and 23. The im- 
portance of upper altitude wind velocity is dramatically demonstrated by Figure 
20 and the accompanying wind profiles (Figure 21) for each day. Wind profile 
for the Washington, D.C. vicinity is measured twice daily by balloon flights and 
is available from the National Weather Record Center, Asheville, North Carolina. 
The wind data presented was taken at 6:15 p.m. on the day of each series of 
stellar observations. On day 269 maximum wind velocity at the tropopause was 
unusually low (approximately ten meters/second), resulting in a sharp reduction 
in scintillation frequency content On the other hand maximum wind velocity on 
day 301 was five times as great. Measured scintillation frequency content showed 
a corresponding increase of a factor of five for this day. The frequency, f a , at 
which normalized spectral density of scintillation is down to one twentieth of its 
zero frequency value can be used as a measure of bandwidth. Data presented in 
Table 1 1 for f B and maximum wind velocity reveal almost a linear dependence 


Table 11 

Scintillation Bandwidth Dependence on Wind Velocity 




Maximum Wind Velocity 

Altitude 

Night j 

f B (Hz) 

(meters/second) 

(kilometers) 

269 

35 

10 

13 

292 

80 

21 

13 

301 

142 

49 

10 


of bandwidth on wind velocity. Scintillation data was taken with a six inch aper 
ture and at approximately the same zenith angle. 

The zenith angle dependence of frequency content is presented in Figure 22. 
Data was taken for two stars with the same aperture: Capella, near the hori- 
zon and other, Deneb, near zenith. Bandwidth, f B , for these two data runs 
taken only nine minutes apart reveal a square root of secant (zenith angle) re- 
lationship predicted by theory. Finally the aperture dependence of frequency 
content is demonstrated in Figure 23. Data was taken for the star Capella with 
apertures of six inches and one inch. Data runs were taken only 21 minutes 
apart thus effectively removing any confusion due to changing zenith angle or 
wind velocity. Bandwidth reduction is a factor of 2 for the larger aperture. 

CIV values measured with the one inch and six inch aperture were .39 and .24 
respectively. This agrees well with Protheroe's measured reduction of a factor 
of 2.5. 

Optical transfer functions were measured for the stars Sirius (zenith angle 
= 70 °) and Capella (zenith angle = 50°) on November 18, J968. These are plotted 
in Figure 24 along with the predicted OTF (Equation 17) using Hufnagel's average 
turbulent profile data (Figure 10). Average image profiles were constructed for 
each OTF by the two-dimensional spatial Fourier transform technique described 
in Appendix A. The results plotted in Figure 25 are characterized by an rms 
image blur circle. This is a measure of image profile radius. Rms values 
exhibit the expected zenith angle dependence. 
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Representative image motion data is presented in Figure 26 for night 292 and 
includes observations over several hours for stars in all parts of the night sky. 
The plot of rms image motion versus secant (zenith angle) includes the theoretical 
square root of zenith angle drawn through the value of image motion at zenith. 

Data agrees well with this dependence, although data scatter increases for large 
zenith angles. Spectral densities of image motion, Figure 27, reveal that es- 
sentially all motion is confined to a bandwidth of 30 Hz. This contrasts with 
bandwidths of 200 Hz or more for irradiance fluctuations. Further reduction is 
apparent in frequency content for large zenith angles compared to zenith view- 
ing. Comparison of image motion and profile data with that obtained by other 
researches is risky at best. However, data of Hosfeld and Johnson et al. can be 
used for order of magnitude comparison. 
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APPENDIX A 

IMAGE INTENSITY PROFILE AND OPTICAL TRANSFER FUNCTION 

The diffraction image formed in an optical system is characterized by a 
normalized intensity distribution called the optical spread function, S. 

S = S(x, y; x\ y'; k) - S(x, x';\) (Ai) 

(light flux in (\ + d\) which originates in dx' and passes through the optics to 
fall on dx, divided by dx dx' d\) 

x' = two dimensional vector in object plane 
x = two dimensional vector in image plane 
The optical spread function is assumed to be time independent and normalized 
such that: 

J" S( x; x ' ; k) dx r 1 

Use of a transfer function to describe the optical system requires the exist- 
ence of an isoplanatism patch, area in the object field over which wave aberra- 
tions are effectively constant. This enables the optical spread function to be 
written: 

S = S(x-x';\) (A 2 ) 
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Equation A4 for S and Equation A7 for I allow image intensity to be expressed 
in terms of the transfer function. 


I(x) = S(x - x' ; X) O(x') dx' 

‘' — CD 

- J*J" T(?> exp [2v i 7 • (x- x' )] df <Xx')dx' 

<x> 

~ U 6XP ‘ (* " *' )] 0(x') dx 1 df 

0(x' ) exp (- 27i i 7 ■ x' ) dx' df 


(AS) 


■r 


T(?) exp (2 tt i f • x) I C 

J-eo 

= T(?)exp( 27 ri f *x)W(?)df 

J -03 


(A9) 


T( f ) is in general complex and can be written in terms of a modulus and com- 
plex argument. 

T(f) = M(f ) exp [277 i <«f )] 


For real symmetric images the transfer function is real. T( f ) ■ = M( f ) 

1(5) = | M( 7) exp ( 2a i 7 • x) d7 (All) 

A transformation to polar coordinates for both x and f allows I to be expressed 
in a form more suitable for computation. 

l(x, y) - f , g) exp [2 tt i ( fx + gy)] df dg 

i - (1) cos 0 x ~ r cos 0 

g - at sin 0 y - r sin d 

n 2” 

M(w, 0) exp [2 t 7 i r at cos (0 - d) 1 at d0 dot 

) 

Symmetry arguments require that I be S independent and M be <t> independent. 


W(f) = J 

r°> 

S(x' )exp (“ 2v \ f • x')dx* = 1 

r® 


or 

I(r) ^ 

1 M((« ) at I exp ( i r at cos 0) d0 dot 

•'0 ■'0 


II 

{* 

T( f ) exp ( 2 tt i f • x) df 

— CD 

(A10) 



1(0 = I" M(aO^J 0 ( 277 ru)dJ 

•'0 


(A 12) 


This is the final form relating image intensity and the transfer function. 
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APPENDIX B 
WAVEFORM ANALYSIS 

Electrical signals generated by the detector package of the SIM are amplitude 
and frequency modulated by atmospheric optical effects. The electronic pre- 
processing system separately determines the effects on signal modulation of 
changes in stellar image intensity, size, and position. A clearer understanding 
of this process can be achieved by modeling the atmospheric effects with single 
frequency sinusoidal components. The chopping of the stellar image is a further 
modulation that must be considered. The voltage, from the photomultiplier tube 
is represented by: 



V„ = voltage magnitude 
V 0 /2 = D.C. level 

V Q = magnitude of chopping frequency component 
- chopping frequency = 3.1 KHz 
C - phase of chopping frequency component 
The actual PMT voltage includes, in addition, frequencies at integral multiples 
of the fundamental, o> c . These are beyond the frequency range of interest and 
will be ignored here. The voltage magnitude, V ()1 is intensity modulated by 
image intensity variations, assumed for this analysis to be at the single 
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frequency, v A . 

o> A 1 200 Hz . 

V 0 = ( 1 - m A COS “* t ) 

m A = depth of intensity modulation 

Variation in image size do not affect the D.C. value of the waveform but cause an 
intensity modulation of the magnitude of the chopping frequency component, 

v o m o/ 2> 

n 0 = (l-m„cosw,t) 

= image size modulation frequency 
m m = depth of size modulation 

Image position changes are reflected in phase, <t > , of the chopped waveform. 

Phase modulation is assumed to be at the single frequency, 

<t> - »V* in,u * 

The PMT output voltage under these models becomes: 

V(t) = (l -m A cosMj t cos o> m cos (fo. t + m^sin t )j 

In the electronic preprocessing system V(t) is split two ways. The first way 
passes through a low pass filter where all signals at or near the chopping fre- 
quency are blocked. 
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The output of the low pass filter, 

!l % \ 

V L< t > \2 ~ ~T cos t ) 

is just the intensity modulation of the image. The second path is through a band 
pass filter where all low frequency components are blocked. 

M 1 ) = ( 1_ m A c °s^ A t)^5 - -3- cosw B t)cos(ai c t+m^sin^t) 

Both filter outputs are fed into a divider section where the operation 

V„(t)V L (t) = V D (t) 

is performed. This step removes intensity modulation and passes the frequency 
components due to size and position changes as sidebands about the chopping 
frequency. The structure of these sidebands can be examined through trigonometric 
identities. 

V n( t ) 5 ("2 ' T COSu 'm l ) cos (“-C * + sir, “<* 1 ) 

- -J cos t ♦ m* sin ^ t ) 

- ~2 cos u- c t cos sin a>^ t j 

m o 

T sin t sin ( sin t ) 
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while the effect of lower frequency image motion is removed by the cancellation 


v D (t) 


m 0 cos t 



K) c " s ( 2k 1 ) 

k = 1 


- m„ sin (o 


t 2 


^ J 2k*l ("><.) Sin [< 2k * 1 >'■'•/ l ] 


= 7 Jo (m*) cos to, t 


+ 7 * cos t <1 ' ‘ 2k “') t ]} 

k=i 


a 

--J ^J 2 k+1 k){c°s[K -(2k ♦ l),)t] - cos[(„ c ♦ (2k ♦ 

k=2 

HI Q 

- T Jo (">*) costo c t 

‘ T J i ( m « ) [ cos {‘“c -“■•*) t ■ cos (" c + “* ) l ] 

■* ~r J 2 ( ) [ COS (“c ' 2u >4>) t + COS (“c + 2u> *) l ] 

- ~2 J 3 ( m <#. ) [ cos (“c * 3u '*) 1 + COS (“c + 3 “*) ‘] 


of Bessel function terms. Thus the peak value is a measure of image size varia- 
tions, represented by 

m Q - ( 1 - rn m cos to ni t ^ 

from D.C, to 200 Hz and image motion greater than 200 Hz. This corresponds 
to a measurement of the short term blur circle of the image. Short term is the 
reciprocal of 200 Hz or 5 milliseconds. 

Image motion at frequencies lower than 200 Hz is measured separately by 
comparison of phases of chopping frequency components from photomultiplier 
tube and reference signal detector. 


J (m # ) is the Bessel function of order i and argument, m$, the depth of modu- 
lation of image position. V B ( t ) is fed into a peak detector stage where the peak 
values of frequency components in the passband ±200 Hz around 3 KHz are de- 
tected. Image motion at frequencies greater than 200 Hz reduces the peak value. 
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APPENDIX C 


SYSTEM COMPONENT SPECIFICATION 


Manufacturers' specification sheets are reproduced below for the following 


15 MAR '67 A 
D055S H 


AMPEX FR-1300 
INSTRUMENTATION RECORDER 


system components. 


PERFORMANCE SPECIFICATIONS 


1. Ampex FR-1300 instrumentation tape recorder. 

2. RCA - 6199 photomultiplier tube. 

3. Venus Scientific photomultiplier high voltage supply. 

4. United Detector Technology PIN-10 photodiode and preamplifier. 

5. Trygon HR40-3B D.C. power supply. 

6. Electro-Craft E-650 motor speed control system. 


General Description 


Deads 


The Ampei FR 1300 it a compact. lightweight portable recorder, avail 
able in 7 or M channel versions. It tenures an integral capstan serve 
system that assures speed accuracy under variable power conditions 
and pcemdes electrical switching over a ranee ot m speeds For ap 
plications tequiruig precisely accurate speeds or where it is desiref 
to record and reproduce on different recorders under unpredictable 
voltage conditions, a control track generato'. demodulator module can 
easily be added to the integral capstan servo system to provide a com 
plete record 'reproduce tape speed servo control system (Speedlock’ 
Transport contro- circuitry is completely intenocier to prevent tape 
damage ■' the recorder is operated carelessly O' accidentally Inter 
changeability ot electronics is possible with nine other Ampes recorders 
using LS 100 electronics Signal compatibility with earlier recorder! 
having 100 Kc direct o> 10 Re FM capability is also provided A tact 
mounted version is available 


Tape Transport 


Tape Speeds: 6C 3C. Is. 7% 3>A 1’/, ips standard All su transpor 
speeds are selected by a single front panel contro All speeds are 
synchronously controlled by a phase locked servo system on cap 
stan drive motor Other hied or variable speeds on special orde> 


Capstan Speed Accuracy: -t005% maximum. long term, whan using 
the internal crystal reference Frequency standard is accurate to 
*001%. long term 

Tap* Speed Devietiea: t 0.25% 

Reels: The FR 1300 tape transport accommodates 10 1 /, inch reels 
Ampei Precision or NAB 

Tap* Specifications. Available in versons for l / ; or I -Inch tap* of 1 
mil or 155 mil Polyester or 1 55 mil Acetate Stated performance 
guaranteed only when using recommended Ampei Instruments 
lion Tape 

Centrals: Illuminated pushbuttons lot Record. Drive Stop Fonvarc 
and Rewind All functions may be remotely controlled Contro- 
circuitry is completely interlocked so that recorder can tie switched 
tram any mode to any other mode without damaging the tap*- 

Fast Wind Time: For 1055 inch reel, with 3600 teet of tape approi. 
mately 3 0 minutes at 60 cycles per second (3 5 mm it 50 cps) 

Start Tima: Time required Irom start command to meet flutter spec 
licetions is as follows 

Speed (ips) Time 

60 8 seconds 

30 6 seconds 

15 and lower speeds 4 seconds 

Slap Time: Mmmum of 1 5 seconds at 60 ips tape speed, shorter 
stop times with loner tape speeds 


Tape Speed Bandpass 

(ips) (cpsi 

60 0 2 to 10.000 

30 0 2 to 5 000 

15 0 2 to 2 500 

7',, 02 to 1,250 

35, 0 2 to 625 

l'/i 0 2 to 312 


Bandpass Flutter 
(cps) (%) 

0 2 to 312 0 20 

0 2 to 312 0 25 

0 2 to 312 0.4 

0 75 0 2 to 312 0 6 

1 2 0 2 to 31? 1.2 

1 5 0 2 to 312 15 


Flettar: Manmum cumulative flutlar (% peak to- peak) 
Flutter 
<%) 

06 
06 
06 


fiap Scatter Trailing edges for record heads (or gap centers tor repte 
duce needs' within a band 100 microinchrs wide (00001 inch) 
Cap Alimeth: All stacks within ? 1 minute ot arc perpendicular to head 
base peate 

Tract DimeauMs: Track width it 0 050 inch, tape track spacing 007C 
center (IRIC Standard i Other heads on special ordci 
Number af Tracks / on inc K 14 o- l inch flRIT Standard) Othr 
heads on special orde* 

lataratack Spaang I 5 rtO 0005 mu . gap 10 gat 


Direct Record' he pronucr System 


frequency Ntapamt: 

Tape Speed 
(IPS' 


S/N* (du> 
Bandpass 

filmed** Unfittere' 


60 

300 cps to 300 Re : 3 0b 

35 

31 

3d 

150 cps lo 150 M -3 dh 

y. 

30 

1' 

100 cps to 75 Ri - 3 ot- 

3? 

27 

? l /» 

50 cps tr 38 Rc : 3 dt, 

30 

22 

3 'a 

50 cps to 19 Rc * 3 df 

79 

72 

I’/a 

50 cps to 10 Rc * 3 db 

2t 

22 


■t« an all Direct anla* aaly I* a ■>«•* 
duvet S/N aH uvutiis *m, h >I)kIm 


Oiracl/FM ayitaw 


"Mmiww) at witpnt W a Wn«pau til'w twvuif II 4k/a 
anawwiwn kayan« Iranti Mela* 


MIS Srgui-te-Nwu Retie: See tabi- 

Harmewc Dtttertien: lass than 2% total ot a 1 Rc signal raoprted 
at 60 ips 

Input lava): 1 0 volt rms nominal id dbv) to produce normal recording 
level operable trom 0 2 to 10 volts rms by adjustment ot input 
potentiometer 

Inptrt F* pedant* Nominal 50 R onms resistive, in parallel with ISC 
pF, unbalanced to ground 

OwDot level: 10 volt rms nominal (0 dtivj. eciota a 600 ah ms v 
graater impedance 

Ovtpart Impedance, less than 50 ohms, unbalanced to ground 


FM Record /Reproduce System 

Frequency Raapmaae: 


Tape 

Frequency Response 

S/N Ratio 

Harmonic 

Speed 

(within 1 0 dbi 

RMS 

Distortion 

60 ipa 

0 to 70000 cps 

46 db 

1.2% 

30 ips 

0 to 10 00O cps 

46 db 

12 

15 ips 

0 to 000 cps 

45 db 

17 

1 u t 'P» 

0 lo 2 500 cps 

45 db 

12 

3>a ips 

0 It. 1 750 cps 

42 db 

15 

1 ’/a ips 

0 to C76 cps 

40 db 

18 


MS StgeeMa-Neiee lobe (it untir wmf, Set tatiti 
Nermemc BnUttien- Swlebii 

DC Drift" Leu then • 0 5% of tub deviation over a four hour panoc 
after warmuo till mmuies test than 2% m 8 hours with lew 
per ituia variations tiatwte i J V. • anc 1?5‘< 

conhnun.' 
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FM Record /Reproduce System icon’t ) 

Ricord/Rtpraducv VelUgi Linearity: 2:0.75% oi lull band, ol a zero 
based straight line 

Input Lewi: Input ol 1 vo It rm ( 0 dbv) to produce -i4Q% deviation, 
operable from 0 5 to 25 volts rms by adjustment of input poten 
hometer 

Input Impedance: Nominal 20.000 ohms resistive, in parallel with 150 
pF, unbalanced to ground 

Output level; 1 0 volt rms (nominal) into 10 K ohms or greater load 
impedance 

Output Impedance: 600 ohms, unbalanced to ground 


Servo Tape Speed l 'antral S\ stem (con'u 

Time Displacement Error 

Tape Speed tiror 

60 ips 100 microseconds 

30 ips 150 microseconds 

15 ips 200 microseconds 

Synchroniiation Time: Aftei 10 minute warmup 
Tape Speed Time 

60 ips 8 seconds 

30 ips 6 seconds 

15 ips 4 seconds 


PDM Record/ Reproduce System 

(IR1G Compatible ) 


System Pulse Characteristics: 
Tape Pulse Duration 

Speed in Microseconds 
(ips) Mm Max. 


60 20 10.000 

30 30 10,000 

15 60 10,000 

7 Vi 120 10.000 

Input level: I volt peak to peak 
10,000 microseconds pulse 
volts peak to- peak. 


Pulse Accuracy m Microseconds 
Over Specified Pulse Widths 
—2 0 (30 to 900 psec duration) 
-20 (45 to 900 «ec duration) 

—3 0 (60 to 2.000 nsec duration) 
—10.0 (120 to 3.000 nsec duration) 

rectangular wave 20 microseconds to 
duration, operable from 0 75 to 20 


tnput Impedance: 20.000 ohms nominal paralleled by 150 pF maxi 
mum. unbalanced to ground. 


Output Impedance: less than 100 ohms, unbalanced to ground 

Output Rise and Fall Time: 2 microseconds maximum (from 10% to 
90% amplitude level) 


Servo Tape Speed Control System 

(Speedlock) 

Frequency Standard Accuracy: ±0.01%. long term 
Spatd Centctiaft Range: i2% ol nominal. 


Power Requirements 

Voltage: 105 to 125 volts, single phase 48 to 62 cps AC (with fre 
quency Standard) 

Power Contamptiin-. less than 325 wafts steady state tor a 14 track 
record/feproduce system Starting surges up to 400 watts may 
be encountered. 

Environment 

Temperature: Operating 5*C to 52*C (40*F to 125‘F) 

Storage/non operating. -30*C to 66*C (-20* to 150*F) 

Altitude Operating 15.000 feet (4600 m> 

Non operating 50.000 feet (15.000 m) 

Relative Humidity: 5 to 95%. non condensing, both operating and non 
operating 

Vibration: Normal handling and transportation only 


Physical Characteristics 

Size: Portable case 24 inches (61 cm) high, by 18 inches (46 cm) wide, 
by 12V> inches (31 cm) deep, for complete 14 track itccrd/repro 
due e system end capstan servo 

Vertical Rack Space Riquired: Rack mounted version Trensport 24 l /i 
inches (62 2 cm), electronics tray S'/, inches Each electronics 
tray accommodates up to 14 record or reproduce modules, plus 3 
auxiliary modules. 

Weight Portable varuon Approximately HO lbs (50 kg) for 7 trick 
system 




RCA-6199 is a 10-stage, head-on, W-diaineter 
type of photomultiplier tube intended for use in acin- 
H nation counters and foe the detection and measurement 
of low-level radiation. 

The spectral response of the 6199, at the 10-per- 
cent points, covers the approximate range from 3200 
angstroms to 6100 angstroms as shown in Fig.J. Maxi- 
mum relative response occurs at about 4400 angstroms , 
and peak cathode quantum efficiency at about 4200 
angstroms. 

SPECTRAL RESPONSE CHARACTERISTICS 
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10-Stoge Haod-On Type 
$-11 Spectral Response 


• Pari Tima Revolution Character! «tl«» - 
Anode Putia Rite Tima: 

2.5 nanaaoeoed* o» 1250 volt* 

Electron Tramit Tima: 

29 nanotecendt at 1250 valtt 

a High Current AaiplJfleaWan - 
1 x 10 6 at 1000 valti 

a Plot Faceplate lor Eoxe In Sclntllloter Ceupllng 


UAt A 

General: 

Spectral Response S'* 1 

Wavelength of Maximum 

Response 4400 ♦ 500 angstroms 

Cathode. Semitransparent Cesium-AnUmony 

Shape RatrC.rcular 

Minimum area . ... 1.2 in 

Minimum diameter 124 in 

Window. . . Lime Glass. Coming® 0080. or etpuvalent 

Shape Piano-Piano 

Index of refraction at 5893 angatroma 1.51 

Dynodea: 

Substrata Nickel 

Secondary-Emitting Surface .... Cesium- Antimony 

Structure Circular-Cag* 

Direct biterelectrode Capacitances <Ap**ox.): 

Anode to dynode No. 10 4 

Anodo to all other electrodes 7 pF 

Maximum Overall Length 4 *^7 » n 

Seated Length 3.88 ♦ 0.19 in 

Maximum Diameter 

Bulb Tl2 

Small-Shell Duodacal 12-Pin. 

(JEDEC No.B12~43>, Non-hygroecopic 

Socket Eby* Part No.9058, or equivalent 

Magnetic Shield Milled P»t No.80802C, 

or equivalent 

Operating Position 

Weight (Approx. 2.2 ox 

Maximum Ratlin*. Abaoluta-Maximu* Values: 

DC Supply Voltage: 

Between anode and cathode .... 1250 max. volte 
Between »ode and dynode 

No. 10 250 max, volte 
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6199. 


Maximum Ratings (Cant'd) 

Between consecutive dynodes . . . 200 max. volts 

Between dynode No. 1 

and cathode 300 max. volts 

Average Anode Current* 0.75 max. mA 

Ambient Temperature* 75 max. °C 


Characteristics Rang* Values (or Equipment Design: 

Under conditions with supply voltage (E) across 
voltage divider providing 1/6 of E between cath- 
ode and dynode No.l;l/12 of E (oreach succeed- 
ing dynode stage; and 1/12 of E between dynode 
No. 10 and anode. 

With E = 1000 volts (Except as noted) 

Min. Typical Max. 

Sensitivity: 

Radiant 9 at 
4400 ang- 
stroms . . - 

Cathode rad- 
iant* 1 at 
4400 ang- 
stroms . . - 

Luminous* . 10 

Cathode lu- 
minous: 

With tungsten 

liiht L 1 S 

source* 3.0xl(T & 4.5 xlO" 5 - A/lm 

With blue 

light 

source* 2.8x10"° - A 

Cathode quan- 
tum efficiency 
at 4200 ang- 
stroms . . - 

Current Ampli- 
fication ... - 

Equivalent Anode- 
Dark-Current 

Input" j ” 

Anode Dark ~ 

Current"* 9 . . — 

Dirk Current to Any 
Electrode Except 
Anode (At 22°C) - 
Equivalent Noise 

Input' 1 _ 

■Anode-Pulse 
RiseTime*. 

Electron Transit 
Time* 


•Made by Coming Glaus Works. Coming, New York. 

^Made by Hugh H. Kby Company, 4701 Germantown Avenue, 
Philadelphia 44, Hn. 

‘Made by Jinnee Millvn Manufacturing Company, 150 Ex- 
change Street, Malden 46. Maaaachuaette. 

^The maximum ruling* in the tabulated data are established 
in accordance with the following definition of the Absolute- 
Maximum Rating System tar rating electron devieea. Ab- 
solute-Maximum ratings arc limiting values of operating 
and environmental conditions applicable to any electron 
device of a specified type aa defined by its published data, 
and should not be exceeded under the worst probable con- 
ditions. 

The device manufacturer chooses these values to provide 
acceptable serviceability of the device, taking no respon- 
sibility for equipment variations, environment variations, 
and the effects of changes in operating conditions due to 
variations in device characteristics. 

The equipment manufacturer should design eo that initially 
and throughout life no Absolute-Maximum value for the in- 
tended service is exceeded with any device under the worst 
probable operating conditions with respect to supply-volt- 
age variation, equipment component variation, equipment 
control adjustment, load variation, signal variation, en- 
vironmental conditions, and variations in device charac- 
teristics. 

•Averaged over any interval of 30 seconds maximum. 

*Tube operation at room temperature or below is recommended. 

•This value is calculated from the typical value for luminous 
sensitivity using a conversion factor of 804 lumens per watt. 

**Thia value is calculated from the typical value for cathode 
luminous sensitivity using a conversion factor of 804 lumens 
per watt. 

* Under the following cooditiooa: The light source is a tung- 
sten-filament lamp having a lime-glass envelope. It is oper- 
ated at a color temperatire of 2870° K and a light input of 
10 microlumena ia used. 

^ Under (be following conditions: The light source is a tung- 
sten-filament lamp having a lime-glass envelope. It is oper- 
ated at a color temperature of 2870° K. The value of light 
flux ia 0.01 lumen and 167 volts are applied between cath- 
ode and all other electrodes connected as anode. 

"Under the following conditions: Light incident on the cath- 
ode ia transmitted through a blue filter (Corning C.S. No. 
&>S8, Glass Code No. 5113 polished to 1/2 stock thickness- 
Manufactured by the Coming Glass Works, Coming, New 
York! from a tungsten-filament lamp operated at a color 
temperature of 2870° K. The value of light flux incident on 
the filter is 0.01 lumen and 167 volts are applied between 
cathode and all other electrodes connected as anode. 

"Measured at a tube temperature of 22° C. Dark current may 
be reduced by use of a refrigerant. 

•Measured with supply voltage (E> adjusted to give a lumi- 
nous sensitivity of 20 amperes per lumen. Dark current is 
measured with no incident light on tube. 


3.6 xlO 4 - A/W 

0.036 - A/W 

45 300 A/lm 


10 - % 
lxIO 8 

2.3xlOT 10 l> 2.5xlCT 9 r 1m 
2.8xir ,3 0 3.1xllT I2 4 W 
4.5 xlCT 9 - A 

7.5 xl(T 7 A 

4xW 12 1.7xl<r n !m 
5x1(T 15 A 2.1x10T 14 4 W 

2.8 x MT 9 - sec 

3.3 xKT 8 - sec 


W VENUS SCIENTIFIC INC. 

PHOTOMULTIPLIER POWER SUPPLY 

MODEL K 15 Reliable • Efficient • Compact • Versatile 



Eout^lOOxEta • Reverse Polarity Protected 

Rippla < 0.1% Peak-Peak • ShortCircuttProtaeted 
Output isolated from Input • All Silicon Components 


The K15 is a high voltage power supply whose versatility 
makes it a very useful component in a number of ap- 
plications. The unit >$ fully encapsulated and shielded in 
a nickel-plated steel can The two output high voltage 
leads are floating and are fully isolated from the two 
input leads. This permits operation from voltage sources 
having either plus or minus polarity while still allowing 
the user to employ the unit as a plus or minus high 
voltage power supply in addition, the case is isolated 
from both input and output and can be grounded by the 
user at a remote point The high voltage filter capacitor 
inside the unit is 006 ul at 3 XV. This keeps the ripple 
at full load down to 0.1%, Where even lest ripple is 
desired, the addition of .01 uf across the output will 
reduce this figure to approximately 0.03%. 

The K15 is protected against improper polarity voltages 
applied lo its input by a series diode. The presence of 
this protective diode accounts lor the most pert for me 
temperature coefficient of output voltage of +0.02%/'C. 
The K1S is 100% tested to guarantee reliable trouble-free 
operation to the user. The electrical oparatmg character- 
istics for a typical unit are shown on the graph for input 
voltages of 3 volts to 15VDC and load* of 10 mega and 
1.25 megs. A plot of Input current as a function of input 
and loads Is also included From this plot it can be teen 
that the operating efficiency of the K15 is better then 
50% for a load of 1 ms. 

Applications: Phot omuttipf left . . . Flying spot scanners 
. . . image dissectors . . . ionization chambers . Scan 
converters . . Geiger tubes . . . Cathode ray tuba bias 
supply . . . Sclnlllation tubes. 

Model K1BU: For ihoae users requiring fast response to 
input voltage changes, the K15 can be supplied without 
internal capacitors this will enable the user to arrive at 
the best compromise between response time and output 
ripple. 


Photo 2 « Actual Sim 


VENUS SCIENTIFIC INC., 25 Bkxxnmidlle Roxd. Hkklvilk, L.I., N.V. 11801 TafcptlOM 51C— 433-33W 
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MODEL K-1 5 

(Electrical Cl»«r»cteri*tlc*) 


c 0 m-vort» oc 



)n -Volt* DC 


physical specifications electrical specifications 


Dimensions 

Volume 

WilfM 

Shock 

Vibration . 

Acceleration 
Thermal Shock 

FWoh, Case 
Can 


1" dia x 2 V4" length 
1 8 cu. m 
3.5 ounce* 

40 g's per MIL-STD-810, Method 516, 
Procedure IV 

20 g s per MIL STD 810. Method 514, 

Curve E, Figure 514-3 

MIL-STD-810, Method 513 

— 54*C to + 71'Cper MIL-STD-810. 

Method 504. Class 2 

Bright Nickel Plate 

cold rolled steel 


K15.A-IUM7 


E output proportion*! to E input (see curve above) 

Floating Output 2500V mtulation strength with respect 
to input ground 

Temperature — 55 X to +7 1'C (operating) 

Temperature coefficient is less then 02%/X 
Efficiency > 50% S full load 

MTBF 105.000 Hr* 7» 71X.(Mil-hdbk 217-A) 

Electrical connections es shown below 





- CKjtKiE | High Vo I Ufi output 

— ClUlt I llkjitln* with itiptcl to Inpull 


TUIOW case 



I O' Sdiccet Cjbbtt lt*tl 


VENUS SCIENTIFIC INC., 25 Btoomingdale Road. Hicksville. L. I . N. Y. 1 1801 T«M»ne 516 — 433-3360 
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FET Model 100 
PREAMPLIFIER 

• Ultra Low Noise 
* Ultra Compact 
• Battery Powered 
Preamplifier 


The UOT Model 100 is an ultra compact, ultra 
tow noise FET preamplifier. The amphtier is 
especially designed tor use with high imped 
ance radiation detectors such as silicon 
photodiodes An internal battery powers the 
amplifier as well as providing alnas voltage tor 
the photodiode It is only necessary to connect 
the photodiode to the input BNC connector of 
the Model 100 amplifier, and an oscilloscope 
or AC vtvm to the amplifier output BNC con 
nector to form an extremely sensitive photo- 
detection system The internal batteries and 
completely shielded case enclosure eliminate 
ground loop hum pick up in the Model 100 
amplifier. Ultra low noise discrete components 
are used throughout the amplifier, yielding a 
wideband input noise voltage of less than 15 
microvolts rms The input impedance is 1 
megohm, giving a wideband noise equivalent 
current input of less than 15 picoamperes, and 
a N E P for silicon photodiodes and narrow 
bandwidth* ol less than 10 12 watts The 
passband is within 3 db from 3 cps to 75 
kc, the gain bandwidth product extends 
beyond 1 me . and the nsetime is 
'■' J y less than 5 microseconds 



UNITED DETECTOR 
TECHNOLOGY 

P.O. Bon 2261 

Santa Monica, California 90405, U.S.A. 
Phona: (213) 393*3785 
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>; jS»S8SB8^^ • ... 

jjj CONNECTION TO PHOTODIODE 

i Input bnc Oscilloscope 





MECHANICAL DETAIL 


CONN: 

BNC Receptacle Jack. 

F 

— 

3.69 

T- 


1st END 

UG290A/U. 




1 88 f 

CONN: 

BNC Racaptada Jack. 





2nd END 

UG290A/U. 

rr— 

T 




Anodlnd 

Finish Die Cast Aluminum 


4 
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ORDERING INFORMATION 

ORDER FROM: UniUd Detector Technology. PO. Box 2251. Santa Monica. California 90405 USA 
SHIPPED: Air Expreat. Prepaid within Continental U S.A PHONE (213i 457 2"i)4 

DELIVERY: 14 daye. TERMS Net 30 days 



FAST, SENSITIVE, 
LARGE AREA 
PHOTODETECTOR FOR 
LIGHT DETECTION 
APPLICATIONS 


PIN-lO’s exclusive combination 
of characteristics, not available in 
other photodetectors, makes it 
ideally suited for a wide variety of 
radiation detection applications. 
Among these are LASER systems, 
spectrometer systems, pattern rec- 
ognition systems, light waveform 
and color detection systems 



Only PiN-10 

offers all these 
features, plus 
reasonable pric- 
ing, each $39. 
ppd. in quantity 
of ten or less. 


The PIN-10 covers three times as many angstrom units as a photo- 
multiplier and has a short wave length response not obtained in 
ordinary silicon pn detectors. Spectral response is from ultra-violet 
to near infrared. 


The PtN-10 is faster than any other solid state detector, as fast as a 
photomultiplier, and 100 times faster than ordinary silicon pn. 
Response time is measured In n»no- seconds. 


Large active area of 3 cm*. 


The PIN -10 offers response from ultraviolet to near infrared rangt, 
yet it has a detectivity as larga as PbS. Tht dark currant is lass than 
a microampare, whllt tha ilght current for a few foot-candies Illumi- 
nation is greatar than a milliampere. allowing the use of PIN -10 
in unchopped light applications. 


Simple circuitry requirements . . . Convenient mounting in one inch 
lens holder . . . BNC connector termination on sealed package . . . 
Order from: 


UNITED DETECTOR TECHNOLOGY 



P.O. Box 2251 
Santa Monica, Calif. 
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E-600 Motor Generator 


Available as a Complete Packaged System or Unwired Component 
System 

MOTOMATIC is a unique* dc servo motor and speed control system designed, engineered and 
manutactured by Electro Craft Corporation. MOTOMATIC is not an SCR type motor speed 
control MOTOMATIC is a complete closed loop system consisting of a fast response perma 
nent magnet motor generator and transistorized velocity feedback amplifier MOTOMATIC 
systems offer precise speed control even in the presence of changing load conditions and 
line voltage 

SPEED RANGE 

MOTOMATIC systems have a nominal speed range of 1000:1. for example, 3 rpm to 3000 rpm 
for a standard E-600 system Compare this with 6:1 for "open-loop" control or 20:1 for better 
quality SCR controls Typical MOTOMATIC speed, torque and output power ranges arc given 
for both E-600 and E-650 systems under SPECIFICATIONS. Standard gearmotor options allow 
an even greater selection of speed and load combinations. Electro-Craft molor-generators are 
quality engineered and manufactured to exacting specifications. Typical brush life is over 
5000 hours at rated speed, as proven by more than 200.000 systems in the field Refer to 
Bulletin No EM-60/65 for motor-generator specifications. 

CONTROLS 

The E-600 M Master control includes speed control, zero to 3000 rpm electronic tachometer, 
torque monitor, adjustable torque limiting, 'eversing-brake switch and remote control capa 
bllities Standard control E-600 S is basically the same as the E 600 M control but without 
tachometer, remote control and torque monitor The E-600 O open chassis version is similar 
to tlie E-600 S except is assembled on an open chassis for installation into customer's 
enclosure The speed control potentiometer, pilot light and switches are supplied loose for 
installation in user's control panel 

Precision E-650 controls are similarly available as master, standard, open chassis and unwired 
systems Heavy duty power supply, high output transistors in the control and high efficiency 
oriented ferrite magnets in the motor, boost the E-650 system capabilities All E-600/650 
series control amplifiers have "»»rn drift" differential Innilt Stan* 


ADDITIONAL Full torque and smooth cogging tree shaft rotation even at lowest 
FEATURES: sp«4s Temperature compensated generator tor long term stability 

Shielded Pall bearings "[QMCv 

W 

nwired Component Uorf 0 . gnd 

Phone 


s drift" differential Input stage 


Numerous specially designed MOTOMATIC systems hBve been manufactured to meet the 
needs of customers whose specifications call for non-standard motor and control configure 
tions Determine your need and let Electro-Craft assume full system design, production, and 
quality control responsibility Please submit complete details of your application to permit our 
engineers to make a proposal 


E-650 MG Motor Generator 



E-650 O Open Chassis Control 


[ SOME MOTOMATIC APPLICATIONS 


Officu Copying Machines 

Cameras 

f'linirigraphii and Prim 

Micro-Film Render-Printers 

Viscosimeters 

Processing Equipment 

Medic. >1 Equipment 

Timing [devices 

Automatic Assembly Machines 

Accurate Motoring Pumps 

Weidmg Feeds 

M.ichuio Tool Programmers 

Laboratory Instruments 

Racer Cnart Drives 

(.tiui[iiiti , i Equipment 

Material Unix llmg Equipment 

Textile Machines 

over Jf»CJ di»f«!».‘Mt appl-Crt* S r i" 
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PERFORMANCE SPECIFICATIONS 

Overall system refutation is represented by the speed torque curvet shown. MOTOMATtC E-600/650 systems 
•re designed for operation within an ambient temperature range not exceeding 0 35 *C end op to 50 *C with 
20% derating. When subjected to load torque changes less than 50% of full load torque a system speed 
regulation for the E 600 of better than 1% may be expected over a large pert of the operating range of 
the system. E-650 systems include a special tow drift amplifier and tow regulation motors. These systems 
maintain an overall speed accuracy of better than 0.5% over the greater part of the system dynamic range 


E-600 SYSTEM PERFORMANCE E-650 SYSTEM PERFORMANCE 







Figure 1. Stellar Image Monitor Experiment Concept 
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Figure 2. Stellor Image Monitor Field Equipment 
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WIDE FIELD 
FINDER SCOPE 


MICROSCOPE 

OBJECTIVE 


SLIDING 

EYEPIECE 


OPTICAL 
HIGH PASS 
FILTER 


STARLIGHT 



6 INCH 
f/8 LENS 


PHOTOMULTIPLIER 
^TUBE 



NARROW FIELD 
EYEPIECE 




REFERENCE 

LIGHT CHOPPER 
DISK 


PHOTODIODE 


CHOPPER DETAIL 


Figure 3. Optical System and Detector Package 
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Figure 5. Electronic Pre-processing System 






STAR 

TAPE LEADER CALIBRATION NUMBER 


DATA 


DEFOCUS BACKGROUND 



Figure 6. Data Run Format 
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Figure 8. Detoil of Electronic Data Collection Equipment 



Figure 9. Typical Data Signals 





Figure 12. Lowposs Filter Frequency Response 


INPUT: AC =1.43 VOLTS RMS 



Figure 13. Bandpass Filter Frequency Respor 



MODULATION INDEX OUTPUT 
(VOLTS) 


A 



Figure 14. Laboratory Calibration of Modulation Index 
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SPATIAL FREQUENCY (CYCLES/ARC SECOND) 
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•3u*pu»d*(] <4jM|*A P U !M < *!* U *Q |w»3»ds •Dur> ! pDjj| JO •|diuo *3 OJ ojoBij 

(ZH) AON3n63dd 



NORMALIZED SPECTRAL DENSITY (Hz 


NORMALIZED SPECTRAL DENSITY (Hz - 



Figure 22. Example of Irradionce Spectral Density Zenith Angle Dependenci 
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Figure 23. Example of Irradiance Spectral Density Aperture Dependence 
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ANGULAR DISTANCE FROM CENTER OF AVERAGE BLUR CIRCLE (ARC SECONDS) 

Figure 25. Image Intensity Profiles 


Figure 27. Example of Image Motion Spectral Density 


NORMALIZED SPECTRAL DENSITY (Hz -1 ) 

pop 



Figure 26. Image Motion Versus Zenith Angle 
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